Influence of Induced Interactions on the Superfluid Transition in Dilute Fermi Gases 



H. Heiselberg 1 , C. J. Pethick 1 , H. Smith 2 , and L. Viverit 2 



O 
O 
O 

(N 

< 

o 



> 
O 

m 

o 
o 
o 



c 

o 
o 



1 NORDITA, Blegdamsvej 17, DK-2100 Copenhagen 0, Denmark 
2 0rsted Laboratory, H. C. 0rsted Institute, Universitetsparken 5, DK-2100 Copenhagen 0, Denmark 

We calculate the effects of induced interactions on the transition temperature to the BCS state 
in dilute Fermi gases. For a pure Fermi system with 2 species having equal densities, the transi- 
tion temperature is suppressed by a factor (4e) 1//3 ~ 2.2, and for v fermion species, the transition 
temperature is increased by a factor (4e)"^ 3_1 ~ 2.2 V ~ 3 . For mixtures of fermions and bosons the 
exchange of boson density fluctuations gives rise to an attractive interaction, and we estimate the 
increase of the transition temperature due to this effect. 
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The study of possible superfluidity in dilute Fermi 
gases has a long history stretching back to the years 
immediately after the development of the BCS theory 
In recent years the subject has received renewed 
attention for a variety of reasons. The first is that super- 
fluidity of nucleons plays an important role in theories 
of neutron stars and of finite nuclei. The magnitudes 
of superfluid gaps are necessary input for calculations of 
transport properties and neutrino emission rates, and for 
modeling the glitches observed in the rotational periods 
of neutron stars. The second is theoretical interest in 
how properties of a Fermi system change as the strength 
of an attractive interaction is varied [H-H. If the sys- 
tem has a two-body bound state, it will behave at low 
densities as a collection of diatomic molecules, while for 
weaker attraction it will behave as a BCS superfluid. A 
third reason is the possibility of observing the transition 
to the BCS paired state in fermion alkali atom vapors 
in traps Q|. Work on this topic has been spurred on by 
the recent success in cooling dilute alkali atom vapors to 
temperatures below the degeneracy temperature H. 

During the 1990's there have been a number of pa- 
pers that calculate the transition temperature, or equiv- 
alently the zero temperature gap, of dilute Fermi systems 
||[7]|^|l(| . The basic physics of these papers amounts to 
summing ladder diagrams with the bare fermion-fermion 
interaction, and expressing the result in terms of the scat- 
tering length for two-body scattering in vacuo. For two 
spin components with equal densities, these papers pre- 
dict a transition temperature T c given by 

fcTc0 = l^-Epe 1 ^^ « 0ME F e^ 2kFa , (1) 
7r e 

for weak coupling (N(0)\U o \ < 1). Here E F = p 2 F /2m F 
is the Fermi energy, pp — hkp is the Fermi momen- 
tum, and in p is the fermion mass. The quantity N(0) — 
mpkp/(2n 2 h 2 ) is the density of states at the Fermi sur- 
face and the matrix element of the effective interaction 
is Uq — 4Trh 2 a/mp, where a is the scattering length for 
two-body scattering, which is negative for an attractive 
interaction. The number 7 is e , where C k 0.577 is 
Euler's constant. 



The above calculation does not take into account the 
effects of the medium on the two-body interaction. Phys- 
ically these processes, which we refer to as induced 
interactions, correspond to one fermion polarizing the 
medium, and a second fermion is then influenced by this 
polarization. This gives rise to an interaction between 
fermions analogous to the phonon-induced attraction re- 
sponsible for pairing of electrons in metallic supercon- 
ductors. The effects of such interactions on superfluid 
transition temperatures have been studied in dense sys- 
tems. In liquid 3 He they are responsible for the ABM 
state being energetically favored close to the superfluid 
transition temperature, whereas in their absence the BW 
state would be the equilibrium one pd[ | . For neutron mat- 
ter calculations predict that they suppress the superfluid 
gap significantly |T^-|ll|] . The corresponding effect in 
a dilute spin-1/2 Fermi gas was considered by Gorkov 
and Melik-Barkhudarov, who found that the transition 
temperature was suppressed by a factor (4c) 1 / 3 ~ 2.2 
compared with the result of Eq. ([!]) [||. In this paper 
we elucidate the physical origin of the suppression, and 
we derive expressions for the transition temperature for 
fermions with a larger spin degeneracy, and for mixtures 
of fermions and bosons. 

The transition temperature is determined by there be- 
ing a solution to the linearized equation for the gap A p , 



£^(p,pO^|^a p ,, 
P ' /? p' 



(2) 



where £ p = p 2 /2mp~ /i is the energy of a fermion relative 
to the chemical potential /i, and = [exp(£/fcT) + 

is the Fermi distribution function. The total inter- 
action is given by 



U(p, p') = £/bare(P, p') + f/ind(P, p')> 



(3) 



where E/bare(p> p') is the bare two-body interaction, and 
£And(p,p') is the induced interaction. Solving the lin- 
earized gap equation is equivalent to finding the temper- 
ature at which the T-matrix obtained by summing lad- 
der diagrams for the repeated interaction of two particles 
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with equal and opposite momenta and with total energy 
2)1 diverges. It is convenient to eliminate the bare interac- 
tion in favor of Tb(p, p'; 2/Lt), the T- matrix for scattering 
in free space of two fermions, each with energy jj,. This 
is given by 

To(p,p';2/i) = C/ b axe(p,p') 

+ ^are(p, P") J~ W, p'; 2fi). (4) 
p" Zt; P" 

Solving Eq. (|4|) for [/bare and inserting the result into (|^) 
one finds 

A p =]TTo(p,p';2 M )^^A p , 
p' ^' 

EAnd(p,P ) — — A r 



p' 



+ ]T To(p,p";2^)^(7 ind (p" I p')^^A p , (5) 
p',p" K v 

Let us first consider the transition temperature in the ab- 
sence of induced interactions. The factor f°(£, P ')/£, P ' falls 
off rapidly for momenta greater than the Fermi momen- 
tum. According to the standard effective range expan- 
sion, the on-shell T-matrix for small p and p' is its value 
for p = p' = and zero energy, To(0,0;0) = Uq, plus 
terms of order p 2 . The latter terms produce higher-order 
contributions than those from the induced interaction, 
and will be neglected here. Thus in Eq. (|J) the gap may 
be put equal to a constant for momenta less than or of or- 
der the Fermi momentum, and the temperature at which 
there is a non-trivial solution to the equation is given by 



l + U 



p 2 dp 1 



27r 2 a 3 £ e«/ feT ^ + 1 



0. 



(6) 



On performing the integration one arrives at Eq. (|l|). 

We now include the effect of induced interactions, 
which we shall assume to be small compared with the 
interaction of two fermions in vacuo. In Eq. (|) the last 
term represents a final-state interaction. The range of 
the induced interaction is of order the spacing between 
fermions or, for the phonon-induced interaction in mix- 
tures, the coherence length. Both these lengths are large 
compared with the magnitude of the fermion-fermion 
scattering length, and therefore in the region where the 
induced interaction is important the wave function of two 
fermions in the presence of the bare interaction, which is 
equal to 1 — a/r, is essentially equal to unity, and final- 
state effects are negligible. One then finds that the tran- 
sition temperature is given to leading order in Ui n d by 



7r e 2 



(7) 



where 



Uind — 



dcosO 



Uind(pFh,p F ti) 



(8) 



is the average of the induced interaction over the Fermi 
surface. The angle between the two momenta is denoted 
by 6>, i. e. n • n' = cos 9. Note that the frequency de- 
pendence of the induced interaction does not enter, only 
its value at the Fermi surface. When the induced inter- 
action becomes comparable to the bare interaction, the 
effects of the frequency dependence of the interaction will 
become important. 




(c) (d) 
FIG. 1. Diagrams for the induced interactions between two 
fermions in different internal states to second order in the 
effective interaction. 



We now apply this result to Fermi systems, and we be- 
gin by considering a Fermi gas with two internal degrees 
of freedom, denoted in Fig. 1 by up and down arrows. 
We shall take the densities of the two components to be 
equal, since this gives the largest gap. The leading con- 
tributions to the induced interaction are represented di- 
agrammatically in Fig. I . One may ask why they should 
have any effect at all in this limit, since they are for- 
mally of higher order in the density than the leading 
term. To understand this, consider the process shown 
in Fig. 1(a), which is the screening of the interaction be- 
tween two fermions by the other fermions. This is of order 
the density of states times the square of the effective two- 
body interaction, and is therefore of order kpaU^. If one 
alters the effective interaction in the expression ([!]) by an 
amount of relative order kp a, it is easy to see that the 
gap is multiplied by a constant factor, as the calculation 
of Gorkov and Melik-Barkhudarov demonstrated explic- 
itly For a contact interaction, the contributions of 
the diagrams in Fig. 1 give 



U in d(pFn,p F h') = U%L(q) 



where 



L,„ = Y ( m(I + ii^!2ln 



1 + 7] 



1-r? 



(9) 



(10) 
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is the static Lindhard function and r\ = q/2kp. Here, 
fiq = p' — p is the momentum transfer in the interaction 
which, as both particles are at the Fermi surface, is re- 
lated to the scattering angle 9 by q 2 — 2k F (l — cos 6). 
Carrying out the integration over angles in Eq. (pf) gives 



U ind = U 2 N(0) 



1 + 2 In 2 



(11) 



and by insertion of this expression in Eq. (R) we find 



T r = 



(4e)i/3 



(12) 



a result implicit in Ref. We therefore arrive at the 
striking conclusion that induced interactions reduce the 
transition temperature by a factor (Ae) 1 / 3 w 2.2 even in 
the low density limit. It is interesting that calculations 
of the superfluid gap in neutron matter at densities at 
which the low-density limit is inapplicable indicate that 
induced interactions reduce the maximum of the gap as 
a function of density by a comparable factor [p^|-^4[ , but 
this is a coincidence. 

The physics of the suppression is best examined by 
expressing the result in terms of the amplitudes for ex- 
change of density and spin fluctuations. In terms of these 
one finds for the interaction in the spin-singlet channel 
for the pair of fermions 



U ind (q) = [U s (q)-3U t (q)}/2, 



(13) 



where U s and Ut are the amplitudes for exchange of spin- 
singlet particle-hole pairs (density fluctuations) and spin- 
triplet ones (spin fluctuations) respectively. The lead- 
ing contributions to the amplitudes are U s — — Ut — 
—UQL(q), and therefore Ui nd (q) — U 2 L(q). One way of 
expressing this is to note that the contributions from the 
diagrams in Fig. l(a)-(c), which are due to exchange of 
spin and density fluctuations with spin projection zero, 
cancel, leaving the contribution from exchange of spin 
fluctuations with spin projection ±1, Fig. 1(d). The im- 
portant process is therefore the suppression of superflu- 
idity by exchange of spin fluctuations, an effect familiar 
in metallic superconductors and liquid 3 Hc p5[ . 

The result above may easily be generalized to v fermion 
species. For simplicity we assume that the two-body in- 
teraction is independent of the species, and is diagonal in 
the species labels. This is true if the interaction is spin- 
independent and depends only on the distance between 
the two fermions. The contribution from the diagram in 
Fig. 1(a) is proportional to v, because all species con- 
tribute to the closed loop, while the other diagrams in 
Fig. 1 are independent of the number of species, since 
there are no fermion loops. Thus the contribution from 
all diagrams is the same as for two species, except that 
the first term is multiplied by v, and their sum is pro- 
portional to v — 3. The transition temperature is given 
by 



T c = (4e)^ 3 - 1 r c0 , 



and therefore for four or more components the transition 
temperature is increased, rather than decreased. 

The results above are not directly applicable to di- 
lute nuclear matter because the nuclear force is strongly 
dependent on isospin: the neutron-neutron scattering 
length is approximately —18.8 fm, while the neutron- 
proton one is 5.4 fm. As a consequence, the induced 
interaction between two neutrons, say, due to excitation 
of neutron particle-hole pairs is of order (18. 8/5. 4) 2 w 12 
times stronger than that due to excitation of proton pairs, 
and therefore the latter contribution may be neglected to 
a first approximation. The fact that induced interactions 
are important for pairing even at very low densities in 
bulk systems indicates the need to investigate how sig- 
nificant these effects are in the outer parts of finite nuclei. 
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FIG. 2. Diagram for the phonon-induced interaction be- 
tween two fermions. The dashed line is the fermion-boson 
interaction, and the wavy line is a phonon in the boson gas. 



We now consider adding bosons of mass mstoa Fermi 
gas with two species. The induced interaction between 
fermions then contains a contribution due to the ex- 
change of boson density fluctuations, analogous to the 
phonon-induced interaction between electrons in metals 
and to the induced interaction between 3 He atoms in di- 
lute solutions of 3 He in 4 He. It is given by (Hjlnj 



u md = Un F x(q,u), 



(15) 



where Ubf — 4irh clbf I^bf (with tubf 
'2 Hi prriB / {niF + tib ) ) is the effective interaction between 
a boson and a fermion, and x(<3S w) is the density-density 
response function for the bosons as a function of the en- 
ergy transfer Tiui. For a dilute Bose gas the response 
function at zero temperature is given by the result of 
using the Bogoliubov approximation, and is 



n B h 2 q 2 /mB 



{huj) 2 -e°(e a q + 2n B UBB)' 



(16) 



where Ubb = 4:nh 2 aBB/mB is the boson-boson effective 
interaction and e° = h 2 q 2 /2rriB- For particles at the 
Fermi surface the energy transfer is zero, and therefore 
the induced interaction relevant for calculating the tran- 
sition temperature is 



find (9,0) 



U 2 

U BF 



1 



Ubb 1 + {hq/2m B sy 



(17) 



(14) 



where s = {tlbUbb / wis) 2 is the sound velocity in the 
boson gas. The average over the Fermi surface of the 
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boson-induced interaction between fermions that enters 
the equation for the gap is 



^H(p F /m BS ) 
Ubb 



(18) 



where H(x) = ln(l + x 2 )/x 2 . Observe that this result is 
independent of the density of fermions in the limit of low 
fcrmion densities. This is in contrast to the induced inter- 
action due to fermion-fermion interactions alone, which is 
proportional to pp. If we again assume that U^J[ is small 
compared with the fermion-fermion interaction, the tran- 
sition temperature is given by Eq. (Q) , where the induced 
interaction includes terms due to fermion-fermion inter- 
actions as well as fermion-boson interactions. In Fig. 3 
this is plotted as a function of the dimcnsionless quantity 
(Ubf/Ubb\Uff\)H(pf /toss). (To make the notation 
uniform we here denote the fermion-fermion interaction, 
which we assume to be attractive, by Uff rather than Uq 
as we did in the earlier part of the Letter.) The param- 
eter U bf /Ubb\Uff\ is a measure of the importance of 
the induced interactions compared with the direct one. 
Since the function H decreases monotonically from its 
maximum value of 1 for x = 0, the largest relative effects 
are achieved for small fermion densities. As Fig. 3 shows, 
the increase in T c due to the presence of the bosons can 
be considerable. 
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FIG. 3. Transition temperature as a function of the pa- 
rameter (U BF /Ubb\Uff\)H(pf /m_gs), for k F \a\ — 0.2 (lower 
curve) and k F \a\ — 0.5 (upper curve). 



Let us now examine a specific example, a mixture of 
the fermionic atoms 6 Li with the bosonic atoms 87 Rb. 
For the boson-boson scattering length we use clbb = 
109an, o-o being the Bohr radius, and for the fermion- 
fermion one a — — 2160ao- The boson-fermion scattering 
length is not known. If we choose the order-of-magnitude 
estimate \<xbf\ = 100a @, U bf /U B b\Uff\ is only 0.18. 
For tib = 10 14 cm' 3 the increase in T c then amounts to 
factors of 2.7 and 1.3 for fcj?|a| = 0.2 and kp\a\ = 0.5, re- 
spectively. If however |osf| were twice as large, these 
factors would change to 17.6 and 2.2. This indicates 



the importance of obtaining information about boson- 
fermion scattering lengths. 

Finally we remark that more detailed calculations are 
required to take into account induced interactions when 
they become comparable in magnitude to the bare inter- 
action. The details of the frequency dependence of the 
induced interactions will then become important. We 
note that corrections to the effective mass of a fermion 
and to the renormalization constant give corrections to 
T c of higher order than those considered here. 
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